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ELECTROMAGNETIC  SIGNALS  FROM  UNDERGROUND  NUCLEAR  EXPLOSIONS 


t>y 

John  Malik,  Robert  Fitzhugh,  and  Fred  Horauth 


ABSTRACT 

Electromagnetic  fields  and  ground 
currents  resulting  from  underground  nuclear 
explosions  have  been  observed  since  the  first 
such  event,  A  few  measurements  have  been 
reported,  but  most  have  not.  There  also  have 
been  some  speculations  as  to  their  origin;  the 
two  most  generally  proposed  are  the  magnetic 
bubble  and  the  seismoelectric  effect.  The 
evidence  seems  to  favor  the  latter  mechanism. 


I .  INTRODUCTION 

Electromagnetic  fields  and  ground  currents  resulting  from  underground 
nuclear  explosions  have  been  observed  from  the  earliest  such  explosion  (1958). 
Despite  over  two  decades  of  sporadic  investigations,  the  data  are  fragmentary 
and  either  incompletely  or  unconvincingly  interpreted.  Some  qualitative 
postulates  of  some  mechanisms  for  generation  of  the  signals^  were  suggested  in 
March  1964.* 

Using  above-ground  electric  field  probes,  loops,  and  current  trans¬ 
formers  on  cables  with  recording  bandwidths  exceeding  30  MHz,  investigators  at 
Sandia  confirmed  the  existence  of  short-duration  signals  of  low  amplitudes 
(<1  V/m  at  1  km),  which  emanated  from  impedance  discontinuities  (equalizers, 

recording  shelters,  cable  splices,  etc.)  rather  than  from  the  well  head. 

2  3 

Among  the  best  late-time  data  are  those  of  Zablocki,  Fitzhugh,  and  Homuth. 
These  data  generally  are  differential  measurements  of  the  voltage  between 
shallowly  buried  electrodes.  Zablocki  used  a  spacing  of  750  m;  Fitzhugh,  6l 


*This  report  is  a  combination  of  two  memos  written  in  1980 — one  assuming  the 
source  to  be  due  to  electric  or  magnetic  dipoles,  the  other  considering  the 
source  to  be  seismoelectric  in  character.  Then  as  now  there  was  no  definite 
conclusion,  although  the  seismoelectric  mechanism  seems  more  likely.  The 
memos  received  limited  distribution  and  received  no  comment. 
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122  ra;  and  Homuth  used  about  150  m.  Zablocki’s  Bilby  data  are  reproduced  in 
Fig.  1;  he  also  quotes  a  scaling  law 

(|e|  =  2.2  •  10^  yV/m)  , 

with  W  in  kilotons  and  R  in  kilometers,  resulting  from  measurements  on  several 
events.  Examples  of  Fitzhugh's  data  are  shown  in  Fig.  2  for  the  case  of  a 
low-yield  event  buried  in  alluvium  and  in  Fig.  3  for  a  high-yield  event 
(Inlet)  buried  in  rhyolite.  (See  Appendix  A.) 

Both  Bilby  and  Inlet  were  fired  in  cased  holes;  the  hole  for  the  low- 
yield  event  was  uncased.  Both  had  canister/rack  support  cables  of  wire  rope 
from  a  strong-back  at  the  surface  and  a  number  of  coaxial  and  multiconductor 
cables  from  the  downhole  device  canister/rack  to  the  surface.  These  slow,  low 
amplitude  fields  appear  to  propagate  by  diffusion  through  the  earth,  not  by 
the  less  dispersive  path  of  to  the  surface,  thence  along  the  surface  to  the 
electrodes.  If  the  source  could  be  approximated  by  elementary  dipoles,  either 
electric  or  magnetic,  both  modes  ought  to  be  observable  with  comparable 
integrals  over  time. 


e  oz  o.«  OjS  om  to  scco^ 

I  I  I  I  t  « 


Fig.  1 .  Record  of  the  horizontal  electric  field 
components  for  the  Bilby  event  at  7.62  km 
south  (magnetic)  of  the  shot  point. 


AxlButhal  Blvctrle  ri«ld,  T/a 


Fig.  2.  Radial  and  azimuthal  fields  from  a  low-yield 
event  emplaced  in  alluvium. 
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Fig.  3. 


Radial  and  azimuthal  fields  from  an  intermediate- 
yield  event  emplaced  in  rhyolite. 
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II.  DIPOLE  MODELS 

Assuming  the  signals  are  due  to  electrical  or  magnetic  effects,  the 

theoretical  base  for  understanding  the  signals  was  formulated  in  1908  by  A. 
4 

Sommerfeld  with  substantial  contributions  from  B.  van  der  Pol  and  K.  A. 
Norton.  Basically,  the  simplified  problem  is  specifying  the  source  moment  and 
solving  Maxwell's  equation  in  a  two-medium  geometry.  The  Navy  (for  underwater 
communications)  and  the  mining  industry  (for  location  of  ore  deposits)  have 

supported  considerable  investigations.  The  most  notable  investigators  are  A. 
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Banos  (University  of  California  at  Los  Angeles),  R.  K.  Moore  and  W.  E.  Blair 

Y 

(University  of  New  Mexico),  and  J.  R.  Wait  (National  Bureau  of  Standards). 
The  excellent  treatise  by  Banos  gives  the  field  description  from  electric  and 
magnetic  dipoles  on  or  in  a  conduction  half-space;  the  summary  relations  are 
incomplete  (as  are  those  of  Moore  and  Blair)  and  contain  only  the  best 
propagation  path  (to  the  surface,  thence  along  the  surface  to  the  observer). 
Other  investigators  that  should  be  noted  are  C.  N.  Vittitoe  (Sandia 
Laboratory)  and  L.  Miller  (Los  Alamos  National  Laboratory),  who  calculated  the 

O 

above-ground  fields  from  an  underground  horizontal  magnetic  dipole. 

Formulation  of  the  integral  equations,  even  with  the  guidance  of  the 
works  of  Sommerfeld  and  van  der  Pol,  has  been  inconsistent  and  fraught  with 
incompleteness  and  error.  The  works  of  Wait  and  Banos  are  very  useful.  In 

5 

particular,  Banos  shows  that  the  horizontal  electric  dipole  is  the  best 
radiator,  with  the  horizontal  magnetic  dipole  as  second  best.  Both  are 
enormously  better  than  vertical  dipoles;  the  electric  dipoles  are  better 
radiators  than  the  magnetic  ones.  Banos  also  shows  that  the  character  of  the 
fields  is  the  same  for  the  horizontal  electric  and  magnetic  dipole  except  for 
the  replacement  of  the  azimuthal  dependence  (j)  by  <t»  +  ir/2,  and  replacement  of 
the  electric  dipole  moment  p  =  I  dl  by  -imk,  where  m  is  the  magnetic  moment 
and  k  =  (imyo)^^^. 

One  of  the  more  informative  developments  of  the  fields  from  a  horizontal 
6-function  electric  dipole  on  the  surface  or  at  a  depth  h  in  a  conducting 
medium  (adaptable  to  a  horizontal  magnetic  dipole)  is  by  J.  R.  Wait  (1960)."^ 
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For  an  observer  on  the  surface  [z  =  0,  but  h  »»  0],  Wait's  Eqs.  (101), 
(102),  and  (103)  reduce  to: 


and 


2  2  2 

The  argument  of  the  Bessel  functions,  I^,  ' ,  is  p  /T  =  (x  +  y  )  ya/8t 
and  T  *=  St/pa.  It  is  useful  to  use  the  asymptotic  expansion  of  the  functions 
for  large  argument.  (Comparison  of  numerical  evaluations  of  the  asymptotic 
expansions  with  the  exact  values  gives  indistinguishable  results,  <1^  error 
for  the  situations  of  interest.)  (See  Appendix  B.) 


I 

o 


(C) 


ii(c) 

2c/2Tr? 


and 


U)  -  21^  (?)  +  (?)  =  21^  -  (2  +  1)1^ 


The  first  term  in  the  braces  of  e  becomes 

X 


/yg  h 

I6v^  p3  t3/2 


2 

The  characteristic  time  is  pah  /H.  This  term  corresponds  to  propagation 
to  the  surface,  then  along  the  surface  with  negligible  further  loss  to  the 
observer . 
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The  second  term  in  the  braces  of  e^  becomes 

[s  -  2 1  ]  /!  f ) 


5/2 


Setting  h  =  0,  this  is  just  one-half  the  expression  for  propagation  from  a 
dipole  in  an  infinite  medium,  and  it  is  thus  the  term  for  propagation  through 
the  medium  along  the  direct  path  frcm  source  to  observer. 

The  asymptotic  expansion  for  e  shows  that  the  only  propagation  path  is 
the  up-and-over  one. 

The  asymptotic  relations  are 


ill  d1 

e  (t)  =  - 

X  ira 


?  P  5/2  , 


e  (t) 

y 


1  3xy 

4 

p  3/2 
(uoh  \ 

e^vrsL 

8p5 

.  /ir  piah^ 

\  ) 

lit  )\ 

e  =  0  , 

z 


2  2  2 
p  =  X  +  y 


2  2.2 
r  =  p  +  h 
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The  expressions  are  written  for  easier  interpretation;  the  integrals  over  time 

of  all  the  square-bracket  terms  are  unity,  and  the  diffusion  times  are  given 
2  2 

by  T  -  par  /4  or  poh  /4,  depending  on  the  propagation  path.  For  a  source  on 
the  surface  (h  •=  0),  the  terms  involving  reduce  to  delta  functions. 

One  way  a  horizontal  electric  dipole  signal  could  be  created  is  through 
an  excitation  of  the  horizontal  cable  runs  or  a  downhole  asymmetry  in  Compton 
and  return  currents  (hard  to  conceive).  The  function  e  has  a 

maximum  at  t  =  2T/n,  and  its  integral  over  time  is  t  r(n/2  -  1). 

For  the  cases  of  interest  here. 


n 

t 

max 

f 

max 

Integral 

3 

2t/3 

0.410 

't/tT 

5 

2t/5 

0.81 1 

T/ir/2 

7 

2t/7 

2.422 

t*3/^/4 

For  a  source  and/or  a  receiver  below  the  interface,  there  are  two 

components  with  comparable  integrals  over  time:  one  is  due  to  propagation 

through  a  distance  (z  +  h);  the  other  is  due  to  propagation  through  an 

infinite  medium  as  noted  earlier.  To  demonstrate  this  separation,  three  cases 

were  computed  using  the  exact  expressions  for  the  Zablocki  conditions  (range  = 

6.7  km).  With  the  source  and  receiver  at  the  interface,  the  fast  component 

(to  this  degree  of  approximation)  would  be  a  delta  function;  to  spread  out 

this  component,  the  source  was  put  at  a  30-m  depth,  with  the  characteristic 
2 

time  as  yah  /4  »  3  ys.  This  is  shown  in  Fig.  4.  Note  that  the  peak  field  is 

-3/2 

near  the  characteristic  time;  the  late-time  fields  fall  as  t  .  To  show  the 
slow  time  component,  the  source  and  receiver  were  located  at  the  interface; 
the  fast  component — a  delta  function  for  these  conditions--is  omitted.  The 
characteristic  time  here  is  related  to  the  receiver  distance;  that  is,  0.14  s. 
If  the  source  were  at  the  depth  of  burial  of  Bilby  (715  m),  the  field 
dependence  for  both  linear  and  logarithmic  time  scales  would  be  as  given  in 
Fig,  5. 
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CV/ft> 


EX  DELTA 


Fig.  Radial  electric  fields  from  a  delta-function  hori¬ 
zontal  electric  dipole  at  30-  and  0-m  depth. 
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The  fields  measured  by  Zablocki  on  Bilby  ought  to  correspond  to  the  e^ 
fields — i.e.,  oriented  in  direction  along  the  cable  runs — but  the  data  do  not 
really  fit  any  of  these  dependences;  the  data  also  show  an  azimuthal  component 
that  ought  to  be  zero.  The  duration  fits  closely  what  would  be  expected  for 
propagation  through  the  medium.  The  peak  is  at  about  10  ms — too  long  for  the 
fast  component  and  too  short  for  the  long  component.  Nor  does  the  shape 
agree — there  is  no  reversal  of  the  fields. 

There  is  also  interest  in  fields  from  a  magnetic  dipole,  notably  in 
evaluating  the  "magnetic  bubble"  model — that  is,  the  expulsion  of  the  static 
geomagnetic  field  by  the  fireball  and  cavity.  As  noted  earlier,  the  fields 
from  horizontal  electric  and  magnetic  dipoles  are  of  the  same  character,  but 
with  (j)  replaced  by  <|)  +  ir/2  and  p  by  -imk  =  -im/iuya. 
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From  Wait's  Eqs.  (96)  and  (97), 

=  L  {p(r^,t)  u(t)} 


and 


P  =  L  {p(r,t)  u(t)} 


and  for  Eq.  (99) 


p(r,t)  -  (-) 


poT 


3/2 


exp  (-  2  rVT)  , 


T  =  8t/pa  , 


p  =  exp  [-  (ups)^'^^  r1 
r 

dP  ,  A/2  „  P 

dF  =  ■  P  "  F  ’ 


11 


or 


~  (yas) 


1/2 


P 


dr  r 


=  L 


dE^u(t) 


dr 


+  ^  L  {p(r,t)  u(t)} 


(-  +  p)  L  {p(r,t)  u(t)}  . 


This  permits  a  formal  solution  for  the  fields  propagated  through  the  medium. 
In  the  asymptotic  approximation,  the  long-term  fields  from  a  horizontal 
magnetic  dipole  of  moment  m  located  at  a  depth  h  below  the  interface  are 


and 


e  =  0  . 

z 

This  is  for  a  delta-function  moment.  For  a  step-function  moment  (u), 
the  fields  are 


e  =  0  , 

xu 


e  »  /  e  ,  dt  , 

yu  7  y5 


and 


e 

zu 


0  . 


The  zeroes 
for  h  =  0. 


?  2 

come  at  t  =  yar  /2  and  yah  /2.  The  relations  are  shown  in  Fig.  6 
(For  most  cases  the  effect  of  the  factor  containing  h  is  small.) 
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X 


_c/p  -1/y 

Fig.  6.  Graphs  of  F(x)  =  (1/2  -  1/x)  x  e  and  J  F(u)  du  as 

o 

functions  of  x. 

The  magnetic  moment  produced  by  the  expulsion  of  a  field  from  a 
sphere  of  radius  R  is 

2ti 

o 

m  =  -  . 

y 

For  Bilby,  the  cavity  radius  was  about  73  m.  Taking  the  horizontal  component 

-■4 

of  the  geomagnetic  field  as  0.23  •  10  tesla, 
m  =  4.5  X  lo"^  Am^  . 

(Note  that  this  assumes  the  conductivity  of  the  melt  and/or  vapor  is 
infinite.)  The  melt  conductivity  is  probably  in  the  range  of  a  few  siemens 
per  meter;  the  vapor  conductivity  may  be  high  and  may  remain  so  for  some 
minutes.  The  thickness  of  the  melt  is  a  few  tenths  of  a  meter;  hence,  the 
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radius  is  very  near  that  of  the  cavity.  See  Appendix  C  for  cavity-radius  and 
rate-of-growth  scaling  relations.  It  seems  probable  that  the  step-function 
relation  is  more  appropriate;  the  time  dependence,  however,  is  not  appropriate 
(Fig.  6). 

7  2 

However,  taking  the  moment  m  as  4.5  *10  Am  [it  was  assumed  as  m6(t)], 

-2 

the  peak  field  under  all  the  assumptions  including  a  =  10  S/m  is 


e 

y 


2  ♦  ^.5  •  ♦  8  ♦  2.08 
ir  y  (7.6  •  10^)^ 
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10  ^  V/m 


Zablocki  observed  2.3  •  10  ^  V/m. 

The  amplitude  and  time  dependence  of  the  expected  fields  from  the 
"magnetic  bubble"  make  it  a  contender  as  the  source  of  part  of  the  observed 
fields  under  some  assumptions  that  may  not  be  viable.  (The  previous  negative 
evaluation  of  these  was  based  on  use  of  formulae  for  vertical  magnetic 
dipole.)  The  observed  fields,  however,  do  not  exhibit  the  proper  radial  and 
azimuthal  dependences.  The  time  dependences  are  those  due  to  propagation,  and 
any  electrical  source  would  exhibit  similar  time  scales.  There  are  also  later 
signals  in  the  recorded  data  that  cannot  be  explained  by  the  "bubble."  These, 
and  perhaps  part  of  the  prompt  signals,  are  discussed  in  Sec.  III.  It  seems 
probable  that  the  signals  are  due  to  more  than  one  mechanism. 


III.  SEISMOELECTRIC  MODEL 

Slow,  low  amplitude  fields  that  appear  to  diffuse  through  the  earth  have 

2  3 

been  observed  by  Zablocki,  Fitzhugh,  Homuth,  and  others.  Zablocki  and 

9 

Keller  proposed  that  the  source  mechanism  is  the  seismic-electric 
effect,  though  the  evidence  was  not  conclusive.  They  also  measured  similar 
signals  on  the  Scooter  event,  an  explosion  of  500  tons  of  chemical  explosives. 
Earlier,  Martner  and  Sparks^ ^  observed  electrical  signals  from  1-  to  20-lb  HE 
charges  detonated  in  bore  holes.  They  attributed  the  signals  to  the 
electroseismic  effect.  Their  results  indicated  that  the  signal  amplitude  is 
proportional  to  the  square  root  of  the  charge  weight.  Extrapolating  their 
data  to  the  nuclear  case,  their  relation  for  the  electric  fields  in  the  earth 
near  the  epicenter  gives  the  relations 

E  =  0.1  V/ra  , 
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where  the  energy  W  is  in  kilotons.  The  amplitudes  measured  by  Fitzhugh  at 
ranges  of  0.5  to  1.0  km  are  consistent  with  the  relation  (Appendix  A). 

The  seismoelectric  effect  was  first  reported  by  Ivanov  (1939),  who 
alleged  that  the  effect  was  due  to  electrofiltration  phenanena  in  moist  soil. 
Frenkel^  ^  presented  a  mathematical  explanation  based  on  the  Helmholtz  theory 
of  filtration  potentials:  "The  particles  of  moist  soil  cire  immersed  in  water 
which  plays  the  role  of  the  dispersive  medium.  The  boundary  surface  between 
these  particles  and  the  water  is  the  seat  of  electric  double  layers,  whose 
aqueous  side  has  a  diffuse  structure.  The  presence  of  such  layers  explains 
the  connection  between  the  flow  of  water  in  the  capillary  spaces  of  the  soil 
and  the  transfer  of  the  surface  electrical  charges  gives  rise  to  an  electric 
field,  in  which  case  electrical  currents  are  compensated  by  the  volume 
conduction  currents.  According  to  the  theory  of  Helmholtz  and  Sraoluchovski , 
the  difference  of  hydrostatic  pressure  Ap  between  two  points  of  the  soil  must 
be  connected  with  a  difference  of  electrical  potential 


AV 


Mirpa 


Ap 


f 


where  ^  is  the  electrokinetic  potential,  i.e.,  the  potential  drop  in  the 
surface  double  layer,  y  the  viscosity  of  water,  and  a  its  electrical 
conductivity.  Since  the  propagation  of  longitudinal  elastic  waves  in  the  soil 
is  accompanied  by  a  variation  of  the  pressure  in  the  direction  of  propagation, 
it  must  also  be  accompanied  by  a  variation  of  the  electrical  potential  which 
constitutes  the  E-effect  discovered  by  Ivanov....  In  the  case  of  Poiseuille 
flow  in  the  liquid  through  a  capillary  tube  with  radius  r,  the  average 
velocity  of  flow  v  is  connected  with  the  pressure  gradient  by  the  relation 


r^  3p  " 

"  8  y  3x  * 

Substituting  into  the  previous  relation. 


r.  2ec  - 

E  s  - 2—  V 

..2 

irr  0 


15 


Using  some  plausible  values  for  the  parameters, 


E  - 


10 


-11 


0.1 


-U  2  -H 
tt(10  10 


-  1  V/m 


9 


where  c  has  been  taken  as  0.1  V.  We  found  no  data  on  the  electrofiltration 
potentials. 

The  result  might  be  expected  to  produce  observable  effects  at  the  main 
shock  front,  at  the  static  water  level  when  the  shock  reaches  that  layer  or 
some  combination  of  reflections  from  the  Paleozoics,  and  at  "shock"  arrival  at 
the  electrodes. 

The  data  obtained  by  Homuth  on  Lowball  (Appendix  F)  seem  the  best  to 
check  the  concept;  those  data  require  more  than  a  dipole  source  model.  A 
sketch  of  geologic  structure  is  given  in  Fig.  7.  These  data  obtained  with 


U9CN 


RESiSTwnnr 


YUCCA  VALLEY 


Fig.  7  Geologic  influence  on  EMP  signals. 
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electrodes  spanning  the  Yucca  fault  at  a  range  of  2.9  km  are  shown  in  Figs.  8 
and  F-1.  Average  sound  velocities  in  the  media  are 
WP  to  surface  1 .8  km/s 

WP  region  2.65  km/s 

Paleozoic  rocks  4.5-6  km/s 

The  depths  of  interest  are 
DOB  =  564  m 
SWL  =  501  m 


P  =  900  m 
z 

Some  times  of  arrival  are 

WP  to  SWL 

WP  to  surface 

WP  to  P  to  SWL 
z 

WP  to  surface  to  SWL 

WP  to  P  through  P  to  fault 
z  z 

WP  to  station 

WP  to  P  through  P  then  to  SWL 
z  z 

Slapdown  at  GZ 


0.024  s 
0.31  3 
0.28  s 
0.59  s* 
0.77  3* 
1 .60  s 
1.19  3* 
1 .62  s 


Times  of  appearance  of  signals  from  Figs.  8  and  F-1  are  0.63,  0.78,  1.0, 
and  1.1  s.  Assuming  electrical  signals  are  generated  at  specific  geologic 
discontinuities  and  then  transmitted  as  electrical  signals  to  the  station, 
possible  regions  appear,  as  noted  with  asterisks.  Better  attention  to 
geologic  structure  and  evaluation  of  shock  arrivals  would  change  the  values 
somewhat  but  there  seem  to  be  some  correlations. 

These  and  other  data  are  given  in  the  Appendices  F  and  G. 
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IV.  CONCLUSIONS 

Considerable  data  exist  for  electrical  signals  emanating  from 
underground  nuclear  explosions;  some,  better  than  others.  A  definitive 
explanation  of  the  late-time  signals  cannot  be  made  at  this  time.  The  more 
probable  source  of  the  signals  is  the  seisraoelectric  effect  proposed  by 

9 

Zablocki  and  Keller.  Data  presented  here  support  that  explanation,  but 
dipole  models  cannot  be  ruled  out  as  contributions.  Data  from  an  explosion  in 
a  well-characterized  medium  with  measurement  of  at  least  radial  and  azimuthal 
components  of  the  fields  at  several  azimuths  and  ranges  could  probably  better 
delineate  the  mechanisms.  Effort  and  equipment  for  such  a  measurement  would 
be  extensive.  The  signals  do  not  seem  useful  as  a  diagnostic  tool  with 
present  information. 
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APPENDIX  A 


EARTH  POTENTIAL,  MAGNETIC  FIELD,  AND  CABLE  CURRENT  DATA 
FROM  A  LOW-YIELD  EVENT  (D--)  IN  ALLUVIUM 


Measurements  herein  were  made  by  R.  Fitzhugh.  The  station  locations  are 
shown  in  Fig.  A-1 . 

Depth  of  burial  326  m 

Diameter  of  hole  48  in.,  uncased 

Number  of  cables  -100 

Earth  conductivity  0.01  to  0.02  S/m 

Earth  potential  plates  were  of  lead,  -2  ft  square,  and  buried  to  a  depth 
of  about  1  m.  Typical  circuits  are  given  in  Fig.  A-2.  Data  are  given  in 
Figs.  A-3  through  A-5. 

Magnetic  fields  were  measured  with  50-turn  coils,  1  m  in  diameter. 
Planes  of  the  coils  were  either  vertical  or  normal  to  the  radius  vector  to  the 
working  point;  that  is,  they  measured  the  vertical  and  radial  magnetic  field. 
These  data  are  given  in  Figs.  A-6  and  A-7. 

Currents  on  cables  are  given  in  Figs.  A-8  and  A-9. 

The  recorder  was  a  multichannel  tape  deck  with  -40-kHz  bandwidth. 
Polarities  may  be  uncertain  in  some  cases. 

A  summary  of  these  data  is  given  in  Table  A-I. 
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TABLE  A-I 

SUMMARY  OF  DATA  FROM  THE  LOW-YIELD  EVENT  D-- 


Orientation 

Location 

Fitzhugh's 

Fig.  No. 

Data 

R(m) 

Field 

Conversion 

D— : 

NS 

w 

A-3 

480 

-0.0164 

EW 

w 

A-3 

480 

+0.0164 

NS 

E 

A-4 

415 

-0.0164 

EW 

E 

A- 4 

415 

-h 

-0.0164 

NS 

S 

A-5 

225 

-0.0164 

EW 

S 

A-5 

225 

~^<l> 

-0.0082 

Normal 

S 

A-6 

225 

I"rI 

255 

Vertical 

S 

A-7 

225 

l%l 

255 

Inlet: 

EW 

N 

6 

1054 

-0.0164 

NS 

N 

3 

1054 

-4 

-0.0164 

Diffusion  Times:  (a  = 

0.01  S/m) 

Station 

Distance  (m) 

2 

T  =  yar  /4  (ms) 

D— : 

GZ 

326 

0.334 

West 

580 

1  .06 

East 

516 

0.84 

South 

375 

0.44 

Inlet: 

GZ 

818 

2.1 

North 

1334 

5.6 

D— : 

The  E,  at  the 
<P 

the  S  location 

W  and  E  locations  ought  to  be  zero,  as  should 

if  MD. 

Inlet:  E_  ought  to  be  zero  at  the  N  location. 

n 


Peak 

+0.MM 
-0.50 
+0.0it 
+0.44 
+0.44 
+0.25 
11  g/s 
8  g/s 


+0.024 

+0.032 


the  Ej,  at 

n 
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I 


Fig,  A-1 .  Station  layout  for  the  low-yield  event  D — . 
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CABLE  CURRENT  CiRCUiTS 


P  TO  TRAILER 


DOWN  HOLE 
CABLE 
FROM  QZ 


m 

m 

TO  5:1  ATTEN 


A-2 


Typical  differential  circuits  for  earth-current  measurements 


250  X  volts 


225  m  S 


Fig.  A-6.  Radial  magnetic  field  from  south  location 


Fig.  A-7.  Azimuthal  magnetic  field  from  south  location 


APPENDIX  B 


USEFUL  LAPLACE  TRANSFORMS  AND 
MODIFIED  BESSEL  FUNCTION  RELATIONS 


I.  LAPLACE  TRANSFORMS' 

f(3) 

-k/s 

e 

r-  -V/S 

✓s  e 

1  -k/s 

—  e 
s 

1  -k/s 

—  e 

/i 

-3/2  -k/s 
s  e 


Ref.  No.* 

#82 

#83 

#84 

#85 


*Transform  number  from  Standard  Mathematical  Tables ,  12th  Ed.  (Chemical 
Rubber  Publishing  Company,  Cleveland,  Ohio,  1956). 
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II.  ASYMPTOTIC  EXPANSIONS  OF  MODIFIED  BESSEL  FUNCTION  RELATIONS  FOR  LARGE  X 


Io(x) 


I^(x) 


Ko(x) 


K^(x) 


/2irx 

X 

e 

/zirx 


('■ 


+  • 


3 

8x 


(8x)‘ 

15 


2!(8x)' 


e'""  ( 1  -  1-  + 
2x  ®  V  8x 


I_  e’*  h  +  ^  - 

2x  ®  I  ^  8x 


2!(8x)‘ 

15 

2!(8x)' 


+  •  ■ 


+  • ' 


Relations  for  derivatives  of  modified  Bessel  functions: 


Iq'(x)  =  I^(x) 
Kq»(x)  =  -  K^(x)  . 
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APPENDIX  C 

SCALING  RELATIONS  FOR  UNDERGROUND  NUCLEAR  EXPLOSIONS 


Cavity  Radii 

Dry  alluvium 

Dry  tuff 

Saturated  tuff 

Rhyolite 

Melt  Radius 


17  m 


1/3 

15  W  ^  m 


(13  to  15)  m 


(9  to  11)  m 


4  m 


Vaporization  Radius 


2  W 


1/3 


m 


CEP  Formula  for  Cavity  Radius 


_  70.2 

c  “  r  .^1/4 


m 


M 


(h  in  meters,  p  in  Mg/m'^) 


Early  Cavity  Radius  Growth 


R 


c 


7  wO-2  t°-5 


m 


(t  in  ms) 


Note:  In  all  the  above  formulae,  W  is  in  kilotons. 

See  Fig.  C-1  for  plots  of  radius  vs  time  for  a  1-kt  explosion  at  675  m  in 
alluvium. 
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RADIUS,  m-kt 


TIME,  ms 


Fig.  C-1. 


Radii  vs  time  for  a  1-kt  explosion  at  675  m  in  alluvium.  Scaled 
from  10-  and  300-kt  calculations  by  T.  Cook  of  Los  Alamos.  Initial 
conditions:  energy  in  1-m-radius  iron  gas  sphere  of  mass  6  Mg, 


„l/3 


r  -  r,  w'^3 
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APPENDIX  D 


GEOLOGIC  SETTING  FOR  EVENT  LOWBALL 


TABLE  D-I 


PHYSICAL  PROPERTY  SUMMARY 


Hole:  U7av 


Physical  Property 


Depth 

(m) 


Event:  LOWBALL 
Data 

Averages 


1 . 

Average  Overburden 
Density 

WP  to  surface 

1 .88  Mg/m^ 

2. 

Average  Density* 

WP  +21.3,-61 .0 

1 .91  Mg/m^ 

3. 

Grain  Density 

WP  +  21.3,-61.0 

2.47  Mg/m^ 

4. 

Average  Seismic 
Velocity 

WP  to  surface 

1807  m/s 

5. 

Average  Seismic 
Velocity 

WP  +0,-6l.0 

2649  m/s 

6. 

Average  Water 

Content 

WP  +21.3,-61.0 

21 .4/22.7$* 

7. 

Porosity 

WP  +21 .3,-61 .0 

39/40$** 

8. 

Saturation 

WP  +21.3,-61 .0 

104/108$** 

9. 

Gas-Filled 

Porosity 

WP  +21 .3,-61 .0 

-1/-3$** 

10. 

Average  CO2 

Content 

WP  +21 .3,-61 .0 

<0.5$ 

11 . 

Depths  and  Per  Cent 
of  Swelling  Clay: 

<10$  montraorillonite 

^H^O  corrected. 

**Water  bias  removed/raw  data. 
tPropagated  error. 
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U7ou  U7av  U7z-I  U7ac 

(BEND  IN  SECTION)  (BEND  IN  SECTION)  (PROJECTED)  (PROJECTED) 


OUJCLI-X 


E 


Fig,  D-1,  E-W  geologic  cross  section  through  U7av, 


N-S  geologic  cross  section  through  U7av. 


APPENDIX  E 


EARTH  MOTION  DATA  OBTAINED  BY  R.  FITZHUGH  FROM  EVENT  LOWBALL 


jELOcmr 


APPENDIX  F 


DATA  OBTAINED  BY  F.  HOMUTH  ON  UNDERGROUND  EVENTS 


seconds  (relative) 

LIKW 


Fig.  F-1.  Lowball  Event  EMP — Buried  plates  across  Yucca  fault. 
Distance:  2.9  km 
Azimuth  from  GZ:  N5°W 
Plate  Separation:  150  ra 
Plate  Orientation:  E-W 

(Tangential) 
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seconds  (relative) 

RNXf 

Fig.  F-2.  Draughts  Event  EMP“Buried  plates  near  seep. 
Distance:  4.7  km 
Azimuth  from  GZ:  N64°E 
Plate  Separation:  30  m 
Plate  Orientation:  N35W 


microvolts/meter 


Fig.  F-3. 


Rummy  Event  EMP — Plates  in  drill  hole  UelK. 


Distance:  4.7  km 
Azimuth  from  GZ:  S7“W 
Place  Separation:  494  ra 
Plate  Orientation:  Vertical 
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microvolts/meter 


seconds  (relative) 


Fig.  F-M.  Draughts  Event  EMP — Plates  in  drill  hole  UelK. 
Distance:  5.2  km 

Azimuth  from  GZ:  S39°W 
Plate  Separation:  555  m 
Plate  Orientation:  Vertic 


Sensors  below  water  table 


Surface 

SWL 

^  - 1. 

1 

1 

1 

\ 

1 

1 

1 

1 

i 

1 

1 

i 

u 

1 

1 
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microvol ts/meter 


first  pulse 


n** 


1(1 


ia.  i«IB  BliBBBBBBBMlLBtiUljWBBBBlIBBBBBBBBI 


BBIf  I  iBBBBBBBBBBBBBBB  RBI 

II'  il 


BBI VBBBBBBBBBBBBBBI 


IWBBffJBBBI 
IIIBB  ^P.BBBBI 

ii.<««.l:fmbbbi 

ir  JBJBBBBI 
IBI.BBBBBBBI 


BBIIBBBBI 
IBIIBBBBI 


nil 


nil 


nil 


t.i 


nwcti 


seconds  (relative) 


Fig.  F-5.  Draughts  Event  EMP — Buried  plates  across  Yucca  fault 
Distance:  km 

Azimuth  from  GZ :  N33°W 

Plate  Separation:  150  m 
Plate  Orientation:  E-W 


II 


Fig.  F-6.  Draughts  Event  EMP — Buried  plates  near  U9cn. 
Distance:  4.7  km 
Azimuth  from  GZ:  N4°E 
Plate  Separation:  300  m 
Plate  Orientation:  N30E 
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APPENDIX  G 

FITZHUGH’S  LOWBALL  DATA 


A  number  of  earth  potential  measurements  were  made  by  R.  Fitzhugh  on 
event  Lowball.  The  locations  of  some  of  his  buried  plates  are  shown  on  the 
map  of  Fig.  G-1 ;  ranges  to  the  plates  from  ground  zero,  from  his  trailer  park, 
and  from  each  other  in  some  cases  are  given  in  Table  G-I.  Measurements  were 
made  of  the  plate  potential  compared  with  the  reference  potential  at  the 
trailer  park;  and,  in  some  cases,  measurements  were  also  made  of  the  potential 
difference  between  selected  plates.  A  selection  of  some  of  those  data  is 
depicted  in  Figs.  G-2 — G-6;  the  upper  waveform  in  each  figure  is  on  a  0.03  s 
time  base;  the  lower,  a  6-3  time  base.  Figures  G-7  and  G-8  are  on  3“  and  30- 
ms  time  bases;  Fig.  G-9  is  on  a  30-ms  time  base.  The  figure  number  below  each 
record  is  Fitzhugh 's  designation  from  his  memo. 

TABLE  G-I 

DISTANCE  TO  POTENTIAL  PLATES 


Distance  (km) 

GM 

-  GM-1 

1.37 

GM-4 

1 .04 

GM-6 

1 .48 

GM-8 

2.01 

GM-9 

2.20 

Trailer 

park  -  GM-1 

0.30 

GM-4 

0.57 

GM-6 

1 .20 

GM-8 

1 .02 

GM-9 

1  .40 

GZ 

1 .08 

GM-6 

-  GM-8 

2.05 

GM-4 

-  GM-9 

1 .94 

49 


Fig.  G-1 .  Locations  of  ground  potential  plates  used  on  Lowball. 
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.KILL1VC4.T$  2.P  ^  .HTLU70LTS 


Fig.  G-8.  Differential  between  GM-6  and  GM-8  (upper)  and 
GM-4  and  GM-9  (lower) — 6-3  time  bases. 


Fig.  G-9.  Potential  on  trailer  park  ground — 30-ms  time  bases 


It  is  not  clear  if  the  data  are  readily  interpreted  as  the  trailer  park 
ground  to  reference  ground  showed  a  large  signal  as  well  (Fig.  G-9).  At  least 
some  parts  of  the  data  have  significance.  A  summary  of  information  relating 
to  Lowball  and  to  two  of  the  stations  is  given  in  Table  G-II.  The  initial 
pulse  may  be  due  to  a  signal  propagated  to  the  surface,  thence  along  the 
surface  to  the  plate.  The  pulse  at  about  5  ms  may  be  due  to  the  signal  from 
diffusion  through  the  medium  to  the  trailer  park;  that  time  is  about  5  ms  and 
it  shows  on  all  records  as  it  is  the  reference  level.  As  the  GM-1  plate  is 
near  the  trailer  park,  with  a  diffusion  time  of  7  ms,  it  probably  has  a  mixed 
signal  of  reduced  amplitude.  The  signals  on  the  shorter  time  scale  records 
thus  appear  to  be  diffusion  signals.  They  could  be  due  either  to  the  magnetic 
bubble  or  to  the  electroseismic  effect.  The  azimuthal  dependence  ought  to 
delineate  which,  but  it  does  not.  The  bubble  signal  ought  to  be  less  affected 
by  the  inhomogeneities  of  the  medium  and  exhibit  a  dipole  azimuthal 
dependence.  As  it  does  not,  it  is  discounted  as  a  plausible  mechanism.  The 
electroseismic  signal  ought  to  be  radial,  but  inhomogeneities  could  readily 
introduce  treinsverse  components. 

The  longer  time  base  records  clearly  show  signals  before  seismic  arrival 
at  the  plates.  One  of  these  signals  seems  to  be  simultaneous  with  the  time  of 
shock  arrival  at  the  surface  (0.3  s).  It  appears  at  a  time  corresponding  to 
propagation  to  the  Paleozoics,  thence  back  to  the  depth  of  the  static  water 
level  (0.5  s).  Plate  GM-1  displays  a  signal  at  about  0.7  s,  which  corresponds 
both  to  seismic  arrival  at  the  trailer  park  and  to  the  time  to  the  P^  layer, 
through  the  layer,  thence  up  to  the  SWL.  This  time  is  0.86  s  for  plate  GM-9, 
which  does  show  a  signal  near  this  time.  This  postulate  seems  to  provide  an 
explanation  for  the  signals  that  are  recorded  at  a  location  before  seismic 
arrival  itself.  The  best  evidence,  however,  is  Homuth's  data  (Appendix  F), 
also  from  Lowball. 
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TABLE  G-II 


LOWBALL  SUMMARY 


Depth  of  burial 

564  m 

Medium 

Saturated  tuff 

Static  water  level 

500  m 

Depth  to  Paleozoics 

900  m 

Average  seismic  velocity  WF 

»  to  surface 

1 .8  km/s 

Seismic  velocity  at  WP 

2.65  km/s 

Seismic  velocity  in  P 

5-6  km/s 

Cavity  radius 

65  ra 

Average  conductivity 

0.01  S/m 

Diffusion  times:  t  = 

ua1^/4 

WP  to  surface 

1.1  ms 

WP  to  GM-1 

7  ms 

WP  to  GM-9 

1 6  ms 

WP  to  trailer  park 

5  ms 

Seismic  arrival  times: 

GZ 

0.33  s 

GM-1 

0.83  s 

GM-9 

1 .26  s 

Trailer  park 

0.68  s 

WP  to  P  +  P  to  plate  sub-point  + 
z  z 

P  to  SWL 
z 

GM-1 

0.7  s 

GM-9 

0.86  s 

WP  to  P  +  P  to  SWL: 
z  z 

0.49  s 

Slapdown  time 

1.62  s 
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